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Abstract. We analyzed the power density spectra of a 
sample of 9 neutron star and 9 black hole binaries in the 
low/hard spectral state. In the power density spectra of 
accreting neutron stars with a weak magnetic field a sig- 
nificant power is contained at frequencies close to one kHz. 
At the same time, most Galactic accreting black holes 
demonstrate a strong decline in the power spectra at the 
frequencies higher than 10-50 Hz. We propose to use this 
empirical fact as a method to distinguish the accreting 
neutron stars from black holes in X-ray transients. The X- 
ray transients that demonstrate significant noise in their 
X-ray flux at frequencies above ^ 500 Hz should be con- 
sidered neutron stars. We propose to explain the observed 
difference as a result of the existence of a radiation dom- 
inated spreading layer on the neutron star surface (Inog- 
amov& Sunyaev 1999). The possible very high frequency 
variabilities of this layer are discussed. 

Key words: Accretion, accretion disks - Instabilities - 
Stars:binaries:general - Stars:fundamental parameters - 
Starsmeutron - X-rays: general ~ X-rays: stars 



1. Introduction 

Among the Galactic X-ray sources black holes (BH) distin- 
guish themselves by the shape of their X-ray spectrum (see 
e.g. Tananbaum et al. 1972, sunyaev & Truemper 1979, 



White, Kaluzicnsky fc Swank 1984 , 

Tanaka fc Shibazaki 199(]| )^ In the low spectral state 
BHs emit a significant part of their luminosity at en- 
ergies of hundreds keV (see e.g. Sunyaev et al. 1991 ), 
while neutron stars (NS) radiate much smaller part 
of their total luminosity in this energy range (see e.g. 
Barret et al. 2000). A soft component that is present in 
the spectra of BH binaries in the high(soft) state has 
a characteristic temperature which can be significantly 
lower than that in the spectra of NS binaries with similar 
luminosity. Besides, the BH binaries in this soft spectral 
state demonstrate a hard power law tail (likely without a 
high energy cutoff up to 500-600 keV, see e.g. Sunyaev 
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et al. 1988, 1992; [Kroeger et al. 1996D , whereas bright 
accreting NSs have never yet shown such spectra (e.g. 
Tanaka & Shibazaki 1996). These spectral properties 
were frequently used as a criterion to determine the 
nature of the compact object. We will not discuss here 
the widely accepted methods of establishing neutron 
star systems through the presence of pulsations or X-ray 
bursts. The detection of coherent pulsations indicates 
the presence of the strong magnetic field and rotation of 
the NS. X-ray bursts (type I) demonstrate that nuclear 
explosions occur in the matter that was collected at the 
surface of the NS during the accretion. However there 
is a significant number of sources, mostly transients, for 
which neither X-ray pulsations nor X-ray bursts have 
been observed. Below we propose another method of 
determining the nature of the compact object based on 
its power density spectrum (PDS) at the high frequencies 
f> 10-100 Hz. 



2. Observations and data analysis 

We used the publicly available data of the Proportional 
Counter Array (PCA) aboard the Rossi X-ray Timing 
Explorer (RXTE) obtained in 1996-1998. We chose the 
observations of X-ray binaries in the low/hard spectral 
state with - 10^^ - lO^'^erg/s ~ 0.01 - O.lLcrit, where 
icrit is the critical Eddington luminosity, made in 1996- 
1998. We only used the data when all 5 Proportional 
Counter Units (PCUs) were operational. The count rates 
for the sources vary from a few hundreds counts/sec/PCA 
to several thousands counts/sec/PCA for both NS and BH 
sources. Our sample of NS and BH binaries is presented 
in Table |l|. 

We constructed the power density spectra of the 
X-ray sources using ^^244 /xsec (2~^^ sec) time resolution 
light curves divided into parts containing 8192 bins. 
The obtained PDSs were normalized to the square of 



fractional variability (e.g. Miyamoto et al. 1991) and the 
ideal Poissonian noise component (Ppoiss(/) — 2/i?tot) 
was subtracted from them. Then the obtained PDSs were 
averaged. Fractional rms normalization is very useful 
for our purposes because it has counting statistics noise 
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Fig. 1. Broad band power spectra of X-ray binaries in the low spectral state. Dashed line on the panel of 
Terzan2/lE1724-3045 shows power spectrum of Cyg X-1 scaled according to their mass ratio (/cygX-i x 7 ^ 
/scaled CygX-i)- This Scaling is not enough to delete the discrepancy in the power spectra of these sources. 
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Table 1. The used observations of NS and BH binaries 



Source 


Proposal 


Dates 


Ref. 


Neutron Stars 


GX 354-0* 


P10073 


Feb. 22- Mar.l, 1996 




/low state/ 








GX 354-0* 


P10073 


Feb. 15, 1996 




/high state/ 








4U0614+091 


P30056 


Mar. 13-20, 1998 


1 


4U1608-522 


P30062 


Mar.31-Apr.2, 1998 




/low state/ 








4U1608-522 


P30062 


Mar. 24-27, 1998 




/high state/ 








SAX J1808.4-3658 


P30411 


Apr., 1998 


2 


lE1724-3045/Ter2* 


P10090 


Nov.5-8, 1996 


3 


GS 1826-24 


P30054 


Feb.- June 1998 




4U 1705-44 


P20073 


Apr.l, 1997 




SLX 1735-269 


P20170 


Feb.- 


4 




P20089 


-Oct. 1997 




KS 1731-260 


P30061 


2-6 Oct. 1998 


5 


Cyg X-2 


P30418 


2-6 July 1998 




Black Holes 


Cyg X-1* 


P30157 


Dec. 1997- Feb. 1998 




GX 339-4 


P20183 


1997 


6 


GS 1354-644 


P20431 


Nov. 19-22 


7 


GRS 1915-^105 


P20402 


Feb.9, 1997 


8 


GRO J1655-40 


P20402 


Aug.l4, 1997 


9 


4U1630-47 


P30172 


May- Jun. 1998 


10 


XTE J 1748-288 




July 1998 


11 


GRS 1758-258 


P10231 


1996- 


12 




P10232 


-1998 






P20166 








P30149 






1E1740. 7-2942 


P10231 


1996- 


12 




P10232 


-1998 






P20166 








P30149 







- the data from these observations were also used for anal 
ysis o f the radiation sp ectra 
[1] -|Ford et al. 1997| , [2 



1999 



[3] - plive et al, 
[5] - IMuno et al. 200' 

IZL 



Wijnandsfc van der Klis 1998 
4] - [Wijnandsfc van der Klis 1999a 
[6] - INowak, Wilms fc Dove 1999 



Revnivtsev ct al. 2000 



PrudoJ^ubov^^^hurazovfc Gilfanov 199E 
' 199C| 



Remillard et al. 



TIOT 



11 



M. 



Tomsick fc Kaaret 2000 



Revnivtse v, Trudolyubov fc Borozdin 2000b 



Main et al. 1999 



T12T 



component, with a deviation from the ideal Poisson level 
(-Ppoiss(/) = 2/i?tot) which is fairly independent of a 
source count rate (see Eq. To subtract the real (not 
ideal) noise component due to the counting statistics from 
our PDSs we used the stationary Poissonian noise com- 
ponent level, modi fied by the deadtimc effects described 
in the papers of Vikhlinin, Churazov fc Gilfanov 1994, 



Zhang et al. 1995| , [Zhang ct al. 199q and successfully 



tested in e.g. Morgan, Greiner fc Remillard 1997, 
Nowak et al. 1999| and Jernigan, Klein fc Arons 2000[ 

The explicit form of the used frequency dependent 
counting statistics noise component, modified by deadtime 
effects appears as follows: 



Pdtif) 



2e 



^tot ^tot \ 

N -1 

-^RphTd 



2i?phrd 1 



cos (27rib/) 



2th J 



sin (ttTvIc/) 



(1) 



here / - the frequency; i?tot - total (in our case - for 5 
PCUs) count rate of the source in the energy band of our 
interest; i?ph - count rate of the source in one PCU in 
energy band of our interest, in our case i?tot ~ 5ii!ph; -Rvie 
- count rate of Very Large Events (VLE) in one PCU; 
Td - deadtime value for any detected event, (e.g. 
Jahoda et al. 1996| , |jahoda et al. l"997| ); T'vic ^ VLE preset 
window, 61/is (level 1) or 150/is (level 2) (see piles 1995 ); 
th - time binning in our light curves (w 244/is); N — 
is the number of bins in the single light curve segment 
{N — 8192 in our case); ^ - is the dimensionless parame- 
ter (very close to 1.0) that was introduced to roughly take 
into account that the deadtime corrections for each PCU 
can be not strictly identical and the model assumption of 
stationary Poissonian process for the observed light curve 
is not exactly valid. But in reality ^ value is very close, 
within a few %, to 1.0. The closeness of ^ to 1.0 can be 
the indication that all used model parameters like td and 
Tvie are close to the real ones. 

Following the conservative approach adopted in 



Jernigan, Klein fc Arons 200C we fitted our PDSs at the 
highest frequencies (/ > 600 Hz for NSs and / > 50 Hz for 
BHs) with the model for the white noise component (Eq. 
[0 without the first term, which was already subtracted) 
and the model for the source itself, leaving the deadtime 
model normalization parameters ^ and Ryie free. We used 
the simplest power law model for the intrinsic variabil- 
ity of the sources. The deduced white noise parameters 
^ and Rvie were always consistent with the anticipated 
ones: ^ « 1.0 (within a few %) and Ryie ~ 100 - 180 
cnts/s/PCU. In the cases when the Rvie (VLE count rate) 
parameter was very poorly constrained (has large statis- 
tical errors) we have it frozen at the values determined 
from Standard_l mode data. This conservative approach 
helps us to subtract most of the PDS components that 
are not real components of the observed target. However, 
if the power spectrum of the observed source has weak 
and very flat noise at frequencies of 500-2000 Hz it might 
be treated as a noise component and will not be detected. 
In general, the used approach for the deadtime modeling 
still has some uncertainties such as the exact value of the 
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deadtime for the photons t^, the stabihty of the VLE win- 
dow Tvie and so on. Therefore we beheve that our PDSs 
still suffer from a small systematic uncertainty of the order 
of fewxlO"^ (a few percent of PAr(/)) in units of squared 
rms. 

Since we were mainly interested in the high fre- 
quency continuum of the PDSs we paid special at- 
tention to all the known timing features of the 
PCA instrument (see e.g. PCA Group Internet Home- 



page ^ttp://lhcawww. gsfc.nasa.gov/docs/xray/xte/pca ) 
Apart from carefully estimating the deadtime modifi- 
cation of the Poissonian level in our PDSs, we ig- 
nored PCA channels 0-7 in all the analyzed data to 
avoid background events variability (PCA Group Inter- 
net Homepage; E. Morgan, privat e communication; see 
also Jernigan, Klein fc Arons 2000 ) that is especially im- 



portant for a weak NS sources like SLX 1735-269 and for 
black holes because of their low intrinsic variability at high 
frequencies. 

To complete our power spectra with the low frequency 
part (~ 10~^ — 1 Hz) we used Standard-l mode data (0.125 
sec time resolution). After the above procedure the power 
density spectra were renormalized to units of squared rel- 
ative rms multiplied by frequency. This helps to better 
visualize the distribution of the power over the Fourier 
frequencies (e.g. Bclloni ct al. 1997 ). 

In Fig. 1^ we present the power density spectra of 18 
sources with subtracted frequency dependent white noise 
level obtained with the described procedure. 

We will present below a comparison of the radiation 
spectra of Cyg X-1 and two NSs (Terzan2/lE1724-3045 
and GX354-0). To get them we used the same RXTE 
data as for the timing analysis (see Table |^). The reduc- 
tion of RXTE/PCA and RXTE/HEXTE data was done 
using standard FTOOLS 4.2 tasks and the methods rec- 
ommended by RXTE Guest Observer Facility. 

3. Results 

In Fig. |l| (upper part of plots) we present the obtained 
power density spectra of nine X-ray bursters in the 
low/hard spectral state. For all of these sources the na- 
ture of the compact object is known - a NS with a weak 
magnetic field. In the lower part of Fig. |l] we present the 
power density spectra for nine BH candidates. For two of 



them (namely Cyg X-1, e.g. Gies & Bolton 198(: and GRO 
J1655-40, e.g. Shahbaz et al. 1999) we have the mass de- 



terminations which indicate that these sources can not be 
NSs. The other sources are believed to be black hole bi- 
naries based on their spectral properties. 

One can see the striking similarity between the PDSs 
of the NS binaries and BH candidates in the low spec- 
tral state up to '^lO Hz, which was already stated in 
the literature (see e.g. Wijnandsfc van der Klis 1999b , 



power spectra (> 10-100 Hz) we notice that even in the ab- 
sence of kHz quasiperiodic oscillations (QPOs) there is a 
dramatic difference between these two samples. The PDSs 
of NS binaries have significant broad noise components at 
the frequencies 500-1000 Hz, sometimes together with well 
known kHz QPOs. The power density spectra of the BH 
binaries in the low spectral state (the typical state for Cyg 
X-1) rapidly decrease at frequencies above 10 Hz and be- 
come by one-two orders of magnitude lower than those of 
the NSs in the same frequency band. The low frequency 
parts of the P DSs probably scale with the com pact object 
mass (see e.g. Wijnandsfc van der Klis 1999b or Fig. |l]). 
However, it is likely that the high frequency part of the 
PDSs of neutron stars demonstrates the additional noise 
component, see Fig. |l|). 

The observed empirical fact described above allows us 
to propose a new method for establishing the nature of the 
compact object in bright transient X-ray sources which 
appear a few times a year in the X-ray sky. We propose to 
classify sources whose PDS extends up to several hundred 
Hz without a strong decline (in coordinates frequency x 
power vs. frequency''^) as NSs with a weak magnetic field 
even if the kHz QPO components are not detected in the 
PDS. 

We understand that the sample of sources used in our 
analysis is not overwhelming. Moreover, the same X-ray 
sources in the high spectral state usually demonstrate 
power spectra, that may significantly differ (see Fig. |2|) 
from each other and from the low state PDS shown in 
Fig.|l|. In other words, the NS systems in the high/soft 
spectral state could demonstrate the absence of signifi- 
cant variability at high frequencies. Therefore the absence 
of high frequency variability can not be the reason to state 
that the source is a black hole candidate. However, we 
propose that any source that does demonstrate significant 
variability (continuum or QPO) at frequencies close to 1 
kHz should be considered NS. Note that for four of the 
nine X-ray bursters represented in Fig. |l|, kHz QPOs have 
not yet been observed. 

4. On the nature of the difference in PDS 

The fate of accreting matter is very different in the case of 
accretion onto a BH and NS with a small magnetic field. 
In the case of an accreting black hole all the observed ra- 
diation forms in an accretion disk, its corona or in an 
advection flow. In the case of a neutron star with a weak 
magnetic field a significant part of the total gravitational 
energy should be released in the boundary layer (e.g. 
^hakurafcSunyaev 1988 , Popham fc Sunyaev 2000) or in 
a layer of spreading matter on the surface of the NS 



( [nogamov&Sunyaev 1999). Approximately 2/3 of the to- 



Psaltis, Belloni fc van der Khs 1999| , [Barret et al. 20001) 
However, when we look at the high frequency part of the 



tal energy releases close to the star surface in the case 
of Schwarzschild geometry, see ^unyaevfcShakura 1986 . 
Only about ^--^1/3 of the total energy is released in the ex- 
tended accretion disk. If the neutron star rotates rapidly 
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Fig. 2. Broad band power spectra of the black hole (Cyg 
X-1) and NS binaries (GX354-0, Cyg X-2 and 4U1608- 
522) in the soft/high spectral state. The difference with 
Fig. 1 is clearly seen for at least 4U1605-522 and Cyg X-2. 
GX 354-0 exhibits high frequency variability even in the 
soft /high state. 



these values change somewhat and the accretion disk con- 
tribution increases (see e.g. ^ibgatuUinfcSunyaev 199^ . 
The surface of the NS is an additional source of soft 
photons for Comptonization. This might be the rea- 
son why the spectra of black holes in the low/hard 
spectral state are generally harder than the spectra of 
NSs (e.g. [SunyaevfcTitarchuk 1989| [Sunyaev et al. 1991 , 
Churazov et al. 1997, see also Fig. |). 



4-.1. The instabilities originated in the accretion disk 

In the standard accretion disk, the region of the main 
energy release is a resonator for acoustic and other insta- 
bilities with maximal frequencies close to uj ^ Cs/H 
Wk, here H - the height of the accretion disk, Cg - is 
the speed of sound, - Keplerian frequency. In the 
case of accretion on the BH with slow rotation the re- 
gion of the main energy release extends from ~ 5Rg to 
~ 27 Rg, here i?g = 2GM/c'^. For the heights H < R 
and Mbh ~ 10 Mq there should be a maximal frequency 
/ ~ 100 Hz. Secular (Lightman & Eardley 1974) and 
thermal (Shakura& Sunyaev 1976) instabilities may de- 
velop in the standard accretion disk but only in the re- 
gion where the radiation pressure exceeds that of mat- 
ter. The characteristic time scales of these instabilities 
are considerably longer (r ^ S /ctaJk, where a ^ 0.01-1 
is the viscosity parameter, see ^hakurafcSunyaev 1973 ) 
than the sound oscillation time scales mentioned above. 



10 



0.1 



Cyg x-1 



GX354-0 



lE1724-30/Terzan 3 



10 100 
Energy, keV 

Fig. 3. The energy spectra of some of the analyzed X-ray 
binaries. The source GX354-0 (4U1728-34) is shown in two 
spectral states - low/hard and high/soft. Even in the hard 
state the spectra of NS are much softer than the spectrum 
of Cyg X-1. 



One can anticipate that the Vclikhov-Chandrasckha r mag- 
netorotational instability ( Balbus fc Hawley 1991 ) could 
also result in the accretion disk variations on the 
time scales exceeding the orb ital time. An analysis 
of the disk turbulence (e.g. ^fowakfc Wagoner 1995 ) 
showed that a strong rollover above ~100 Hz should 
present in the power spectra of an accretion disk. 
The phenomenological shot noise model (see e.g. 



Terrel 1972, Lochner, Swank & Szymkowiak 1991) and 
the model of coronal energy release variations (e.g. 



Poutanen & Fabian 1999) consider the similar absence of 



strong variability at frequencies higher than 100 Hz. The 



advcction dominated accretion flow (e.g. [chimaru 1977 
Narayan & Yi 1994) are not subject to the thermal insta- 



bility mentioned above. In spite of the absence of detailed 
studies of different instabilities for the advection flows it 
seems quite reasonable to assume that a limit of the order 
of H/cs may be valid in this case as well. Rapid BH ro- 
tation can lead to increased frequency values, but in this 
case the rotation parameter a should be pretty close to 
the critical value 1. If a < 0.5 the influence of the black 
hole rotation is small. 

All instabilities existing in the accretion disk mod- 
ulate the flow of matter onto the NS surface. There- 
fore, we could expect that the majority of the types of 
variabilities we observe in accreting BHs must manifest 
themselves in accreting NSs with characteristic times pro- 



Shakura&Sunyaev 1976, 


Wijnandsfe van dcr Khs 1999b 


Inogamov&Sunyaev 200C 


. However, in the BH systems 
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Fig. 4. The comparison of the power spectra of black 
hole (Cyg X-1) and neutron star (Terzan 2). The dashed 
line shows the power spectrum of Cyg X-1 scaled accord- 
ing to the mass ratio with Terzan 2 (/cygX-i x 7 — > 
/scaled CygX-i)- This simplc scaling is important but insuf- 
ficient to explain fully the difference in the high frequency 
variability of Terzan 2 and Cyg X-1. Note, that slopes of 
power density spectra of Terzan 2 and Cyg X-1 in the high 
and low frequency limits are similar, however the power 
spectrum of Terzan 2 is sufficiently broader that that of 
Cyg X-1. 



the innermost regions of the accretion disk can generate 
the variability at its Keplerian frequencies (close to ~100- 
200 Hz), but their contribution to the observed PDS can 
be small because of the decrease of the energy release close 
to the last stable orbit (e.g. ghakurafcSunyaev 1973 ). In 
the NS case, where the boundary layer can possess all 
accretion disk frequencies, this variability could be more 
pronounced. 

The simplest assumption is that the characteristics fre- 
quencies in the power spectra of the sources scale as M~^. 
This scaling law is valid for e.g. Keplerian frequency in the 
vicinity of the last stable orbit, thermal and secular insta- 
bihties of the accretion disk in the region of main energy 
release, Balbus-Hawley instability. However, this assump- 
tion is not enough to account for the observed difference 
in the high frequency variability of the NSs and the BHs. 
To demonstrate this we compared the PDS of the neutron 
star system Terzan2/lE1724-3045 {Mi - 1.4M0) with 
M2/M1 times scaled PDS of black hole system Cyg X- 
1 (M2 ~ lOM©) (Fig. |l| and |). It is obvious that such a 
scaling is important but insufficient to explain the high fre- 
quency variability of Tcrzan2. It is most likely that there 
is an additional component coming from the neutron star 
surface or boundary layer in its vicinity. 



4-2. The spreading layer and the highest possible 
frequencies in the PDS of NSs 

The inner boundaries of accretion disks around BHs and 
neutron stars lie at similar radii measured in the values of 
the gravitational radii of the central object. However, we 
expect a small energy release from the region with the ra- 
dius smaller than the inner boundary of the accretion disk 
around a BH. In the case of NS accretion there is a solid 
surface and the matter should finally join this surface. We 
know this because of the existence of X-ray bursters. In ac- 
cretion disks matter rotates with Keplerian velocity. The 
velocity of the surface of the NS is 6-3 times lower (see e.g. 
^trohmayer 1999 or van der Klis 200C for reviews). There- 
fore matter must decelerate and release its kinetic energy 
in the narrow spreading layer on the surface of the NS. 
In the case of very low luminosities, less than 0.01 £crit 
(where Lcrit if the critical Eddington luminosity) all de- 
celeration occurs in the narrow boundary layer between 
the innermost region of the accretion disk and the surface 
of the NS. We could consider this layer an atmosphere 
in which accreting matter is losing its kinetic energy and 
radiating it away (see e.g. ShakurafcSunyaev 1988 ). At 
luminosities above 0.01 Lent the radiation pressure does 
not permit such a simple picture. Matter rotates around 
a NS and slowly spirals up the meridian toward higher 
latitudes. In this flow the centrifugal force and the radi- 
ation pressure force balance gravitation with very high 
precision. Effective gravity (the difference between the 
gravitational force, centrifugal force and radiation pres- 
sure force) becomes hundreds of times smaller than the 
gravitational force. The turbulent friction with underly- 
ing denser layers rotating with the surface velocity of the 
star leads to the deceleration of the flow and energy release 
which is radiated away. This approach predicts the exis- 
tence of two broad bright belts equidistant from the star 
equator. These belts are fed partially with the advection 
of t hermal energy from the r egions closer to the equa- 
tor ( InogamovfcSunyaev 1999 ). The luminosity of these 
bright belts exceeds the luminosity of the accretion disk. 

The spreading layer may be considered a thin broad 
box. The flow along the latitude has Keplerian velocity in 
the vicinity of the equator and slowly decelerates to higher 
latitudes. This velocity is close to a hundred thousand 
kilometers per second. The sound velocity Cg within the 
radiation dominated box is of the order of 20 000 km/sec. 
The highest sound frequencies inside the box are of the 
order of h/cg and might reach even 40 kHz in the case 
of low luminosities and 6-7 kHz in the case of luminos- 
ity close to the Lcrit (Inogamov&Sunyaev 2000), here h is 



the thickness of the spreading layer. Simple consideration 
shows that we could expect a lot of turbulence generated 
sound, resonant sound waves and plasma instabilities with 
characteristic frequencies that do not end at the Keplerian 
frequencies, but can reach several kilohertz. The sound 
waves of much lower frequencies might also exist in the 
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bright belts: one can estimate the characteristic frequen- 
cies / ~ Cs/(27ri?) ~ 300 Hz for the waves propagating 
along the latitude and / ~ Cs/{R/S.9) ~a few kHz (where 

- is the angular width of the bright belts) for the waves 
propagating along the meridian within the belts. 

As well as we expect the Keplerian frequency on the 
last stable orbit (~ 200Hz for Mbh ~ IOM0) to be the 
maximal possible frequency for the BH variability in X- 
rays, 40 kHz and 7 kHz mentioned above might be the the 
maximal frequencies for the NSs PDS. In both cases (for 
BHs and NSs) the power spectra should drop sufficiently 
towards these frequencies, but might be high below them. 

In the spreading layer we have very strong shear in the 
flows along the longitude. Under such circumstances the 
long living whirls like those of thyphoons, Jupiter red spot 
or "cats eyes" due to Kelvin-Helmholtz instability might 
regularly appear. The evolution of the chaotic magnetic 
field structure inflowing into the spreading layer through 
the neck of the accretion disk (region with (Kl jdr = 0) may 
produce bright spots in the spreading layer. The accretion 
disk screens from us half of the star, making it completely 
unobservable. Part of the star surface is not visible due to 
eclipses by the star itself at any inclination angle except i 
close to 0. The rotation of the bright spots leads to their 
regular eclipses. The duration of these eclipses is close to 
half of the period of the matter rotation in the spreading 
layer. The quasi regular eclipses (the velocity of the bright 
spot can slowly change) can also give rise to the additional 
variability (noise) of X-rays from the NS at frequencies of 
hundreds of Hz. 

There are many different mechanisms of instabilities 
in the narrow (with heights from 400 m up to 1.5 km de- 
pending on the accretion rate or luminosity) , thin spread- 
ing layer. Unfortunately, at the present level of sensitivity 
of RXTE we can not see a significant variability of the 
sources at the frequencies above 2-4 kHz. It is likely, that 
future missions with higher collecting area of the detec- 
tors can shed light on the very high frequency (/ > 5 — 10 
kHz) variability phenomena. 
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